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Ultramagnetized neutron stars or magnetars are magnetically powered neutron stars. Their strong magnetic 
fields dominate the physical processes in their crusts and their surroundings. The past few years have seen 
several advances in our theoretical and observational understanding of these objects. In spite of a surfeit of 
observations, their spectra are still poorly understood. I will discuss the emission from strongly magnetized 
condensed matter surfaces of neutron stars, recent advances in our expectations of the surface composition of 
magnetars and a model for the non-thermal emission from these objects. 
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1. INTRODUCTION 

Simply put magnetars are neutron stars whose mag- 
netic fields dominate their emission, evolution and 
manifestations. In the late 1970s and early 1980s, a 
fleet of sensitive detectors of high-energy radiation un- 
covered two new phenomena, the soft-gamma repeater 
and the anomalous x-ray pulsar. Strongly magnetized 
neutron stars provide the most compelling model for 
both types of object, and observations over the past 
few years indicate that these phenomena are two man- 
ifestations of the same type of object. Soft-gamma 
repeaters exhibit quiescent emission similar to that of 
anomalous x-ray pulsars [e.gQ, 0, 13, 0, and anoma- 
lous x-ray pulsars sometimes burst What makes 
magnetars a hot topic of research is the rich variety 
of physical phenomena that strong magnetic fields ex- 
hibit. 

This article will focus on the quiescent emission 
from these interesting objects rather than the bursts 
themselves (The reader may wish to refer to the sem- 
inal work of Thompson and Duncan for details of 
the burst but may also want to look at Heyl and Hern- 
quist 1^ for an alternative). Furthermore, the article 
will concentrate on recent results. 



2. MAGNETAR MANIFESTIONS 
2.1. Soft-Gamma Repeaters 

The first soft-gamma repeater was discovered on 
March 5, 1979 when gamma ray detectors on nine 
spacecraft across our solar system recorded an intense 
radiation spike The burst of gamma rays origi- 
nated from near a supernova remnant known as N49 in 
the Large Magellanic Cloud. The tail of the gamma- 
ray burst exhibited an eight-second pulsation (in con- 
trast with the classical gamma-ray bursts which show 
no periodicities). If one combined this eight-second 
period with the age of the supernova remnant and as- 
sumed that the object (presumably a neutron star) 
was born spinning much faster, one estimated a mag- 
netic field on the object of ^ 10^*^ G, much larger than 
any neutron star discovered up to that point. 



This first SGR became known as SGR 0526-66. 
During 1979, two others were discovered: SGR 1806- 
20 (27 December 2004 event), SGR 1900-1-14 (28 Au- 
gust 1998 event) [3- In 1998 SGR 1627-41 became 
the fourth SGR to be discovered |Tll |. 

Thompson & Duncan 0, ^2 argued that the evolu- 
tion of a ultrastrong magnetic field could explain both 
the outbursts and quiescent emission from soft-gamma 
repeaters, and the term "magnetar" was born. They 
argued earlier that if protoneutron star was born spin- 
ning sufficiently rapidly, a dynamo could dramatically 
amplify the standard pulsar magnetic field (~ 10"'^^ G) 
that the protoneutron star was born with to ~ 10^^ G 



or more 



2.2. Anomalous X-ray Pulsars 

IE 2259-1-586 was the first anomalous x-ray pulsar 
to be discovered 0|. In the early nineties these ob- 
jects began to form a unique class 0, 0, 0, ITsj . 
They were found to have much in common with the 
soft-gamma repeaters. 

These objects typically have pulsed X-ray emission 
with steadily increasing periods of several seconds, X- 
ray luminosities of ^ 10'^'^ — 10"^® ergs s~^, soft spec- 
tra, and no detected companions or accretion disks. 
Furthermore, they are typically observed through hy- 
drogen column densities of ^ 10^^ cm~^, indicating 
that they are not common. 

Two of the five confirmed AXPs are located near 
the centres of supernova remnants IE 2259-1-586 
and IE 1841-04514, 181. as well as the AXP can- 
didate AX J1845-0258 :i9!|. The remaining objects 
are 4U 0142-^61 20^, IE 1048.1-5937 21], IRXS 
J170849.0-400910 22] and possibly XTE J1810-197 


The emission from the AXPs fits neatly within the 
magnetar model T^. Heyl & Hernquist „24] argued 
that the thermal flux passing through an ultramagne- 
tized, hydrogen or helium envelope is sufficient to ac- 
count for the x-ray emission from these objects. Hey 
& Kulkarni psj examined how magnetic field decay 
can augment the thermal energy budget. For fields 
less than about 10^^ G magnetic field decay in a real- 
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istic model does not strongly affect the emission from 
young (less than 10,000 years) AXPs but can greatly 
increase their lifetime as observable x-ray sources. 

Alternative models such as accretion have fallen by 
the wayside because even very low mass companions 
have not been discovered orbiting these neutron stars 
nor has the tell-tale optical emission from even a trun- 
cated accretion disk been detected. 

2.3. Strongly Mangetized Radio Pulsars 

Moreover, a number of radio pulsars have been dis- 
covered with inferred magnetic field strengths similar 
to those of magnetars and apparently exceeding the 
value Bqed ~ 4.4 x 10^^ G 26, 27, 28] It is not clear 
why these objects have magnetic fields comparable to 
the AXPs but do not exhibit AXP-like emission. The 
inferred magnetic field of AXP 1E2259+586 is actu- 
ally smaller than PSR J1847-0130, an otherwise ordi- 
nary radio pulsar. The AXPs and SGRs exhibit much 
greater x-ray emission than their relatively inactive 
cousins. § 13 speculates further on the possible differ- 
ences between these two types of objects. 



3. PHYSICS IN A STRONG MAGNETIC 
FIELD 

Isolated neutron stars generally drawn on several 
sources of energy. For most observed neutron stars, 
the dominant source of energy is the rotation of the 
star. 

E,ot = Im"" « 2 X 1046/45 (-^) erg (1) 
2 \ 1 sec / 

Magnetars are typically young, highly magnetized 
neutron stars. The magnetic energy, 

M = ^ J B^dV ^2x 10^'RlBlss5 ergs, (2) 

exceeds the rotational energy by an order of magni- 
tude. Nearly as large is the thermal energy of the 
star 

U ^ IQ^' RlMl'^Tl^ ergs, (3) 

where Afi.4 and Tg.s are the mass and core temper- 
ature of the star in units of 1.4M0 and 10^'^ K, re- 
spectively, and i?6 is the radius in units of 10 km and 
Sjgg 15 surface magnetic field of the star in units 

of 10^^ G. Much more weakly magnetized objects such 
as AXP 1E2259-I-586 act like magnetars as well. For 
this object U - 1007W. 

It is unclear why some strongly magnetized neutron 
stars behave like "magnetars" rather than just hap- 
pily spinning down and radiating as a radio pulsar. 



Not only do the soft-gamma repeaters and anoma- 
lous x-ray pulsars share a penchant for bursting 
and Kaspi in these proceedings] that the pulsars lack 
they also generally have much stronger x-ray emis- 
sion that radio pulsars ~- they are not only strongly 
magnetized but hot as well. As we have seen, for 
AXP 1E2259-I-586 the thermal energy exceeds a naive 
estimate of the magnetic energy by two orders of mag- 
nitude. It is quite natural to speculate that the heat 
flux travelling through the surface layers of the AXPs 
and SGRs may play an important role in their "mag- 
netar" behaviour. A thermomagnetic interplay [23 
between the strong magnetic field and the strong heat 
flux may help to power their bursts . 

Gamma rays from by far the largest burst from a 
soft-gamma repeater arrived at Earth ten short days 
after the close of this meeting. The recent superbursts 
from SGR 1806-20 [23, El IHl brings many of these 
issues to the fore. Particularly the energ y of this burst 
has been estimated to be ~ 2 x 10^^ erg [33, one hun- 
dred times that of the previous two superbursts from 
soft-gamma repeaters. It has been argued that this 
is a once in a century event |30l ]. Given that detec- 
tors sensitive to such phenomena have only existed 
for 48 years (if one generously assumes that Sputnik 
and the Explorer satellites would have noticed such 
a burst), a once per century rate is rather conserva- 
tive. To account for the number of observed SGRs, the 
soft-gamma repeaters must have a lifetime of several 
thousand years, so they would be expected to exhibit 
several dozens of such bursts. The budget of magnetic 
energy alone is insufficient to account for these bursts, 
so one is forced to conclude that either we are seeing 
the very end of the SGR phase for SGR 1806-20 or 
that the energy reservoir of SGR 1806-20 exceeds the 
magnetic energy stored in the dipole field. 

Subsequent observations of SGR 1806-20 have re- 
duced the estimate for its distance 33], but they still 
only obtain a lower limit to the distance. This result 
may reduce the energy budget for the event and ease 
the energy crisis. Further work is clearly needed. 

Until the past year the observed quiescent emission 
from SGRs and AXPs was thought to be dominated 
by thermal emission. The following subsections dis- 
cuss the various physical processes that determine the 
thermal emission from AXPs and SGRs (and possibly 
the non-thermal emission as well). 

3.1. Vacuum Physics 

For many purposes, vacuum polarization effects can 
be calculated using an effective Lagrangian for the 
electromagnetic field. Following the usual convention, 
we write 

C = + Ci + ■■■ . (4) 

Here, C is the full Lagrangian density, £0 is the classi- 
cal term, and Ci includes vacuum corrections to one- 
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loop order. Higher order radiative corrections would 
be described by additional terms. The second-order 
term is smaller that £i ^34] . 

Because the Lagrangian is Lorentz invariant, both 
Cq and Ci can be expressed in terms of the Lorentz 
invariants of the electromagnetic field 



/ = ^^^,F^'' = 2(|B|' 



lEl 



and 



(5) 



and 



^QED 



B 



Bqeb 



(13) 



where the final equality holds if the electric field can 
be neglected. 

The function Xq can be evaluated analytically (as 
can the higher order terms in the expansion for Ci) 
with the result IstII 



(6) 



where e^''^'' is the completely antisymmetric Levi- 
Civita tensor. The effective Lagrangian of the elec- 
tromagnetic field was derived by Heisenberg and Eu- 
ler |35j| and Weisskopf [HI using electron-hole theory. 
In rationalized Gaussian units, we can write Cq and 
£i as 



(7) 



OO 



dx 

y3 



X 

2 coa{J+ x) + cos( J- x) 



r2 

-°QED 



cos( J+ x) ~ cos( J_ x) 



(8) 
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where 



J± = 



25, 



QED 



-I ± i- 



1/2 



(9) 



a = /he is the fine structure constant, -Bqed = 
m?c^ /eh K, A A x 10^"^ G, and a similar quantity can 
be defined for the electric field, -Eqed = m?c^ /eh k, 
2.2 X 10^^ V/cm. The above expressions for £o and Ci 
are identical to the cor resp onding terms in eq. (45a) 
of Heisenberg & Euler [33 ■ 

The above integral cannot be evaluated explicitly, 
in general. Heyl & Hernquist [33 have derived an 
analytic expression for Ci as a power series in K: 
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(14) 



where 



C = 21n47r - 41nyl --ln2 2.911785285, (15) 
3 

and the constant In A is related to the first derivative 
of the Riemann zeta function, (^^\x), by 

InA = - C^^'(-l) = 0.248754477. (16) 

The integral of Inr(x) can be expressed in terms of 
special functions (eqs. 18, 19 in |33)- (See also Dit- 
trich et al. 1979; Ivanov 1992, but note the cautionary 
remark in \3T\.) 

The expression above for Xq (x) can be expanded in 
either a Taylor series in the weak field limit, ^ <C 1, or 
an asymptotic series in the strong field limit ^ ^ 1, 
as can the higher order terms in equation H1Q|I , to give 
series expansions for Ci as a function of either / and 
K, or equivalently B and E. In particular, to lowest 
order in the weak field limit (^ <C 1) 



1 



a 



90n Bl^^ 



[{B^ - E^f + 7(E-B)2] 



(17) 

In the limit of an ultrastrong magnetic field, B ^ 
Sqed, but for a weak electric field, E <C Eqeb, 
can be written 



£i = £1(1,0) + K 



dCi 
'dK 



K=0 



where the first term in this series is 
'1' 



A(/,0) = ^IXo^^ 



a 
47r 



(10) 



(11) 



-u/i^(.UCOthu+l+^ 

If' V 3 



(12) 



OTT 



In 



— ) 

Bqed J 



12 InA + ln2 



(18) 

(see e.g. eq. 29 in |37| for the higher order terms). We 
note that equations (|17|) and (|18|l agree, respectively, 
with the correspondin g te rms in eqs. (43) and (44) of 
[ssj . The analysis of |37j generalizes the expressions 
of Heisenberg & Euler to arbitrary order. 

Also, note that while our expression for Cq is iden- 
tical to eq. (4-120) of Itzykson & Zuber [^^i equation 
((T7|) diflters from their eq. (4-125) by a factor of l/47r. 
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as a consequence of a difference in tlie system of units 
employed.^ Our expressions for Cq and Ci both dif- 
fer from those in Berestetskii et al. [s^ by an overaU 
factor of l/47r (their eqs. 129.2 and 129.21); however, 
the dynamics of the fields are invariant with respect 
to rescalings of the Lagrangian. 

We emphasize that the expression for the La- 
grangian in the weak- field limit, equation (|17|l . can- 
not be applied to magnetar fields which are thought 
to have i?NS ^ Bq^d. The use of the weak- field ex- 
pressions to calculate, e.g. the index of refraction of 
the vacuum near the surface of a magnetar will result 
in estimates that are incorrect by more than an order 
of magnitude at the relevant field strengths. In this 
limit, the Lagrangian should instead be approximated 
by e.g. equation H18|l . which is an asymptotic series 
for Ci valid for B » i?QED and E <C Eqed- 

After integrating out the effects of the virtual 
electron-positron pairs, one is left with an effective La- 
grangian and many electromagnetic phenomena can 
be treated by simply treating the vacuum itself as a 
medium endowed whose properties are described by 
the Lagrangian. This approximation is valid for pho- 
ton energies much less than TOgC^. For higher energy 
photons the assumption of a homogeneous field im- 
plicit in the above construction of the effective La- 
grangian breaks down and one must resort to other 
methods [for example liof . 



3.2. Atomic Physics 

The atomic physics of materials on the surfaces of 
magnetars is dominated by the magnetic field rather 
that the atomic nucleus. A useful figure of merit is 
the ratio of the cyclotron energy of an electron to the 
binding electron of an electron to a nucleus of atomic 
number Z 



E, 



cyclo 



i^Rydborg Z"^ jlop-meC^ a'^ Z^ BqeD 

B 



1.1^2 X 109 G 



(19) 



(20) 



where Bqed is defined in H3.ll We see for magnetars 
this ratio is large for Z up to one hundred or more. 
The atomic physics is in this regime has been probed 
extensively for hydrogen and helium [for example l4lL 
|43^ but approximately for larger atoms such as iron 

What remains uncertain is how these atoms interact 
and what are the properties of strongly magnetized 
material in the bulk at low pressure. Is it a solid. 



^Itzykson & Zuber l38il use Heaviside's units in defining the 
Coulomb force; and are smaller in this system than in 
ours by a factor of Att. 



liquid or gas? Is it metallic? These questions are cru- 
cial to understand the details of the thermal emission 
from neutron stars. For further details on the physics 
of atoms and molecules in strong magnetic fields the 
reader is urged to consult the excellent review of Lai 
0. 



3.3. Thermal Physics 

Through much of the envelope of a magnetar the 
electron gas is effectively one-dimensional because 
only the first Landau level is typically filled 45, |4^ ■ 
This restriction dramatically affects the thermody- 
namic properties of the material and allows approx- 
imate yet reasonably accurate analytic treatment of 
the thermal conduction through the envelope of a 
magnetar psj . 

The key parameters are the chemical potential of 
the electrons, the temperature and the strength of the 
magnetic field. It is helpful to define the dimensionless 
quantities, 

C^^,x^ = C^-l,r=^ (21) 

in addition to ^ = _B /Bqed from § 13.11 The number 
density of electrons in the degenerate limit (r ^ C~l) 
is given by 



1 



"e,0 



for a weakly magnetized gas {x ^ 2^) and 



^e, magnetar 



(22) 



(23) 



for a strongly magnetized gas. The symbol A denotes 
the Compton wavelength of the electron, h/{mec). 

The pressure of the magnetized gas also differs from 
the unmagnetized case. The number density of elec- 
trons in the degenerate limit (r <C C ~ 1) is given 
by 



Pe,0 



1 



m^c 



for a weakly magnetized gas (x 3> 2^) and 



e,ma£rnetar 



4^^3- 



[xc - Hx + c)] 



(24) 



(25) 



for a strongly magnetized gas. Although the mag- 
netic field changes the pressure dramatically, the elec- 
tric pressure remains isotropic due to the changing 
mangetization as the electrons are compressed |47| . 

An important comparison is at what density do the 
electron become degenerate in the two cases. Magne- 
tar envelopes are generally cool, t <^ lsoC^l~ 
a;2/2 « T at the onset of degeneracy. The ratio of 
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the densities of the gas at the onset of degeneracy in 
unmagnetized and the ultramagnetized case is given 

by 

>te,magnctar _ ^ — ^ 1 (Of{\ 

ne,o ~2x2"4r>>^- ^^^^ 

The ratio of the pressures at the onset of degeneracy is 
given by a similar factor. For a magnetar, ^ ~ 10 — 100 
and r ^ 10~^, yielding a factor of 100—1000 increase 
in the density, column density and pressure at the 
onset of degeneracy. Magnetar envelopes only become 
degenerate densities larger than 

Pnd/d = 3.92 X 10^ g cm-3^(ln5Typicai)'/^ 

.6/7 -5/7 4/7 -1/7 , . 

^ ^56 ^26 -'off,65s,14 • l^'J 

where ^Typical is a typical value of ^/t in the enve- 
lope, ^Typical ~ 6 X 10'^^. The unmagnetized envelope 
become degenerate at densities ~ 10^ — 10^ g cm~'^ 



4. THERMAL EMISSION FROM 
MAGNETARS 

The thermal emission from magnetars is powered 
by a combination of the residual heat from the super- 
nova explosion (Eq.Oand possibly the magnetic field 
(Eq. El . The accretion alternative predicts more opti- 
cal emission than is observed from these objects [but 
see |49', for an alternative point of view] . If the ther- 
mal energy originates from within the star, the energy 
must first travel through the insulating envelope of the 
neutron star. The conductivity of this layer throttles 
the heat flow and determines the total thermal emis- 
sion from these objects. 

In the non-degenerate regime, the conductivity is 
nearly constant and determined by the outgoing heat 
flux both for mag netars and weakly magnetized neu- 
tron stars [45I lia j. The photons in the extraordinary 
polarization mode dominate the heat conduction, and 
the conductivity itself is strongly dependent on the 
magnetic field. In the degenerate regime the conduc- 
tivity increases dramatically. Along the field lines, the 
strong magnetic field strengthens the heat conduction, 
while across the field lines, the field inhibits it. For 
a fixed internal temperature of the star, the emerging 
flux is approximately proportional to Z^'^-^i?" '* cos^ tp 
where tp is the angle between the local magnetic field 
(B) and the normal and Z is the atomic number of the 
nuclei that comprise the envelope 0i ■ Potekhin 
et al. examined the dependence of the cooling of 
neutron stars (including the superfluidity of their in- 
teriors) on the presence of a light-element envelope in 
further detail, verifying the earlier conclusions. 



4.1. Magnetized Atmospheres 

The study of mangetized neutron stars is simplified 
by the fact that neutron stars atmospheres are ex- 
ceedingly thin. Unfortunately, this is an inadequate 
consolation for the complications that the strong mag- 
netic fields introduce. 

• As discussed in § 13.21 the binding energies of 
even light species such as hydrogen and helium 
may exceed several hundred electron volts in the 
magnetar fields, so magnetar atmospheres are 
probably only partially ionized [see Figs. 1 and 2 
of|^. The light element atmospheres of more 
weakly magnetized stars are typically fully ion- 
ized, but even in weak fields heavier species such 
as carbon, oxygen and others will not be fully 
ionized. The presence of bound species affects 
not only the opacity of the plasma but also its 
equation of state and electromagnetic polariza- 
tion. 

• The opacities for the two polarization modes dif- 
fers by several orders of magnitude, so the po- 
larized radiative transfer must be calculated. 

• Vacuum polarization fS l3.1|l complicates the be- 
haviour of the propagation modes of the plasma 
itself. As radiation propagates through the vac- 
uum resonance where the vacuum and plasma 
contribute equally to the index of refraction, 
the mode that below the resonance had the 
higher opacity now has the lower opacity and 
vice versa. This effect complicates the radiative 
transfer for for radiation sli ghtl y above the peak 
of the thermal distribution ^^2iE3- 

• The opacites depend strongly on the direction 
of that the radiation propagates relative to the 
magnetic field. Where the field is normal to the 
surface only the altitude is important (like in un- 
magnetized transfer) but the distribution must 
be sampled more finely to handle the sensitive 
angular dependence. Where the field is makes 
an angle with respect to the normal (the general 
case), the intensity must be calculated both in 
altitude and azimuth. 

The recent work of Potekhin et al.||5j| may be con- 
sidered an example of the state of the art in mag- 
netized atmospheres. Fig. depicts the results for 
calculations that fully include the effects of partial 
ionization and vacuum polarization as well as spec- 
trum using less realistic assumptions. The shape of 
the spectrum changes dramatically as the calculations 
become more sophisticated. Even from these com- 
prehensive simulations, several important effects are 
absent. First, the magnetic field is assumed to point- 
ing in the normal direction. Relaxing this constraint 
makes the calculations more onerous but similar in 
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_ B=10i3 G 




E (keV) 

Figure 1: The total emergent spectrum from a strongly 
magnetized neutron star atmosphere from the 
calculations of Potekhin et al. 54] . A blackbody 
spectrum is compared with several atmospheric 
calculations of increasing realism. The fully ionized 
model assumes that neutral species are absent. The 
hybrid model includes partial ionization to calculate the 
equation of state and the opacity but assumes that 
polarization modes are those of a fully ionized plasma. 
The partially ionized model includes the effects of partial 
ionization on the opacity, equation of state and 
polarization. All of the models include vacuum 
polarization 



principle and is required to simulate the emission from 
the entire surface of the neutron star. Second, the 
magnetic field is still two orders of magnitude less than 
that of magnetars. As the magnetic field increases the 
difference in the opacities between the two modes in- 
creases and also the material in the outer layers may 
undergo a phase transition Finally the surface 

layers may not consist of hydrogen but rather helium 
or heavier elements (see § 14.4(1 . 

A tantalizing alternative is that magnetar atmo- 
spheres aren't atmospheres at all, but rather the sur- 
face of a magnetar is condensed. The condensation 
temperature may be as high as 10® K for iron in fields 
stronger than 10^^ G and even for hydrogen in fields 
of a few xlO^^ G. In this situation detailed balances 
allows the emergent spectrum from the star to be de- 
termined by understanding how light refiects off of it. 
Van Adelsberg et al. |^ calculate emergent thermal 
spectra from the condensed surface of a magnetar. If 
the surface is smooth, the spectrum deviates from a 
blackbody by at most a factor of two and exhibits only 
mild absorption features that correspond to the ion 
cyclotron rcsosonace and the plasma frequency within 
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Figure 2: The total emergent spectrum from a strongly 
magnetized condensed neutron star atmosphere from the 
calculations of van Adelsberg et al. |5^. 

the surface. The emission for rough surfaces resembles 
that of a blackbody more closely. 

4.2. Magnetospheric Propagation 

Regardless of the detailed composition of the at- 
mosphere of magnetars, the vast difference between 
the opacities in the two polarization modes nearly en- 
sures that the radiation will emerge from the atmo- 
sphere highly polarized perpendicular to the direction 
of the local magnetic field. Pavlov and Zavlin [5d| ar- 
gued that because the direction of the magnetic field 
varies across the surface of the neutron stars, the net 
polarization from the entire visible surface is greatly 
diminished. This argument neglects that the region 
surrounding the magnetar is optically active. Specif- 
ically at the x-ray energies, the the intense magnetic 
field decouples the propagation modes in the magneto- 
sphere through vacuum polarization fji l3.1|) . Heyl and 
Shaviv ^] argued that polarized radiation propagates 
adiabatically through the magnetosphere. The adia- 
batic condition holds within the polarization-limiting 
radius [see ,58. for the plasma analogue] , 




where r is the distance from the center of the star, /i is 
the magnetic dipole moment of the neutron star, and 
(5 is the angle between the dipole axis and the line 
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Figure 3: The total emergent flux from the visible 
surface of the star. The magnetic pole makes an angle of 
60° with the line of sight. The radius of the star is 12 km 
and its mass is 1.4 Mq. The effective temperature at the 
magnetic pole is 10®'^ K. The calculation is for 
/i = lO'''^ G cm^, corresponding to a surface field of 
~ 10^'' G respectively. The solid curve traces the total 
flux. The short dashed curve traces the flux polarized 
perpendicular to the projection of the magnetic moment 
in the plane of the sky. The long-dashed curve traces the 
flux polarized parallel to the projected magnetic moment. 
The heavy curves trace the results including vacuum 
polarization and the light curves neglect it. 



of sight. The observed polarization reflects not the 
direction of the magnetic field at the surface of the star 
but at the polarization-limiting radius; consequently, 
if this radius is much larger than the radius of the 
star, the polarized radiation from disparate regions of 
the stellar surfaces adds conherently an a large net 
polarization results. 

Fig. 131 depicts the polarized spectrum summed over 
the entire surface of the neutron star. The local emer- 
gent radiation spectrum is determined from an atmo- 
spheric model that assumes that the plasma is fully 
ionized hydrogen and neglects vacuum polarization 
within the atmosphere. The bold curves trace the 
results from a calculation that includes vacuum polar- 
ization within the magnetosphere and the light curves 
neglect vacuum polarization. If vacuum polarization 
within the magnetosphere is included the radiation is 
more than 99% polarized near the peak of the spec- 
trum. Without it the total emission is more modestly 
polarized by about 10%. 

In the X-ray regime, the radiation from the entire 
surface of the neutron star is highly polarized for sur- 
face magnetic fields greater than 10^^ G or so. At 



lower photon energies, especially in the optical and in- 
frared, the observed polarization is a sensitive probe 
of the magnetic field and plasma density within the 
magnetosphere as well as the radius of the neutron 
star |59| . Polarimetry of magnetars and neutron stars 
in general will provide a unique new probe of their 
environments and interiors. 

4.3. Magnetic Field Decay 

Thompson and Duncan (l^ argued that the decay 
of the strong magnetic field of a magnetar could ac- 
count for the quiescent emission from the surface and 
examine the gradual evolution of the magnetic field of 
the star. Goldreich and Reisenegger 60] examine sev- 
eral modes of magnetic field decay: Ohmic diffusion, 
Hall drift and ambipolar diffusion. These processes 
have the following timescales: 

iohmic ~ 2x 10"-| (-^) yr (29) 

J 8 V Pnuc / 

ilmbip - 3 X 10^^ yr (30) 

-°12 

''ambip rpQ r,2 ^ ''ambip \'^^J 

8 12 

iHaii - 5 X 10«^p (^) yr (32) 

where L5 is a characteristic length scale of the flux 
loops through the outer core in units of 10^ cm, 
is the core temperature in units of 10^ K and B12 
is the field strength in units of 10^^ G. Ohmic de- 
cay dominates in weakly magnetized neutron stars 
(B < 10^^ G), fields of intermediate strength decay 
{B - 10^2 - 1013 G) via Hall drift, and intense fields 
{B > 10^^ G) are mostly strongly affected by ambipo- 
lar diffusion. Thompson and Duncan |l2l | examined 
the possibility of obtaining an equilibrium between 
neutrino cooling and heating through magnetic field 
decay for B - lO^^ G and T - 10^ K. 

In the case of a magnetar older than a few hun- 
dred years, ambipolar diffusion dominates the other 
porcesses at least in the core of the neutron star. 
Ambipolar diffusion involves motion of the charged 
species (electrons, protons and possibly muons) rela- 
tive the the neutrons. In the motion of the change 
species results a change in the density of these parti- 
cles (if the flow is irrotational) , the chemical equilib- 
rium must be reestablished through weak interactions. 
Otherwise the flow is solenoidal and weak reactions 
are not required. 

Using the results of Goldreich and Reisenegger |63] , 
Heyl and Kulkarni examined how magnetic field 
decay would affect the thermal emission from the sur- 
face of a magnetar and its cooling, including a realistic 
treatment of the strongly magnetized envelope of the 
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neutron star. Fig.0]shows the results of their analysis. 
If the magnetic field is sufficiently strong {B ^ 10^^ G) 
field decay has a dramatic effect on the emission from 
young magnetars, less than ten thousand years old. 
Even without field decay the magnetic field dramat- 
ically increases the heat flux through the surface of 
the neutron star. For intermediate fields {B ~ 10^^ G) 
field decay does not strongly affect the emission of the 
young magnetars, but for fields similar to and stronger 
than 10^^ G field decay increases the lifetime of hot 
magnetars by a factor of several. This dramatically in- 
creases the number of hot magnetars that you would 
expect to see in the Galaxy and possibly the evolution 
of soft-gamma repeaters and anomalous x-ray pulsars. 



4.4. Diffusion 

A crucial input for our understadning of the emis- 
sion from thes surface of neutron stars is the compo- 
sition of their outer layers. As discussed in S I3.3I and 
14. li the presence of light elements in these layers can 
have a dramaticially effect on the total thermal emis- 
siokn from the star as well as its spectrum. The strong 
gravitiational acceleration at the surface of the neu- 
tron star {g ^ 10^"' cm s^-^) ensures that the lighter 
elements will float to the top of the atmosphere; con- 
sequently, one would expect the emergent spectrum 
and flux to be determined by the lightest elements 
present in the outermost fluid layers of the star. If 
any material fell back onto the neutron star during or 
after the supernova, it spallates as it hits the surface, 
resulting in hydrogen and helium nuclei. 

The thickness of a pure hydrogen layer is limited 
the stabilty of protons against decaying into neu- 
trons. If the chemical potential of the electrons ex- 
ceeds mass difference between a proton and a neu- 
tron {( ~ 2.53,2: w 2.33) it becomes energetically 
favourable for protons to combine with electrons and 
form neutrons. These neutrons will quickly bind to 
protons, forming deuterons and ultimately helium nu- 
clei. 

According to equations and [531 this occurs at a 
density of about 1.2 x lO^g cm^"^ for a weakly magne- 
tized gas and about 3.4^ x lO^g cm^"^ for a strongly 
magnetized gas in which f ^ 1. 

The thickness of the hydrogen layer that can be 
supported on a magnetar increases with the strength 
of the magnetic field. The thickness of a helium layer 
is assumed to be limited by the shell flash instability 
that results in Type-I X-ray bursts on neutron stars. 
This criterion yields a maximum column density of a 
helium layer of approximately 10*~^ g cm~^ [e.g.l6l|. 
It is unclear whether the strong magnetic field affects 
nuclear burning of helium, but it does affect hydrogen 
burning |63 |. If a neutron star accretes more than 
about 10^^ g of material after its birth, one would 



expect find moving from the outside inward hydrogen, 
helium and then carbon or other heavier elements. 

These results neglect the fact that although the lay- 
ers are stabily stratified, a thermally excited tail of 
protons will penetrate the layer of helium and possi- 
bly be captured onto the carbon and other nuclei lying 
below. Chang, Arras and Bildsten |Q examined this 
problem in detail. As protons diffuse through the he- 
lium layer, their density drops as a power law in the 
non-degenerate regime and exponetially in the degen- 
erate regime. In a magnetar the electrons become 
degenerate at a much higher column density than in 
a neutron star with a smaller magnetic field (Eq. 127(1 . 
Although diffusive nuclear burning is not terribly im- 
portant over the timescale of a few thousand years 
for weakly magnetized neutron stars, in magnetars a 
thin atmospheric hydrogen layer is consumed as pro- 
tons diffuse through a maximally thick helium layer 
in only a few years. 

The efficiency of diffusive nuclear burning on a mag- 
netar indicates that the atmospheres and enevlopes of 
magnetars likely consist of helium or heavier elements 
rather than hydrogen. Because ionized hydrogen at- 
mospheres has generally be the focus of research into 
magnetar atmosphere (see § 14.1(1 clearly more work is 
needed. 



5. NON-THERMAL EMISSION FROM 
MAGNETARS 

Until the past year the quiescent emission from 
magnetars was thought to be mostly due to thermal 
emission from the surface of the star. Their soft x-ray 
spectra are notoriously well described by blackbody 
models (possibly with a power-law component at high 
energies) [for example |6j|. However, recent results 
from INTEGRAL indicate that the thermal emission 
is just the tip of the iceberg jB^E^E^, and Hurley in 
these proceedings]. The optical emission from AXPs 
is also clearly non-thermal (or at least not from the 
surface) |6^ . 

Thompson and Beloborodov [g^ present two mod- 
els to account for the non-thermal hard x-ray emis- 
sion from the AXP IE 1841-045. The first involves 
the heating of a surface layer of the neutron star to 
kT ~ 100 keV by currents driven in the magneto- 
sphere by the twisting of the magnetic field in the 
crust. The second model focusses on a region far from 
the star (at ten stellar radii) where the electron cy- 
clotron resonance is approximately 1 keV. The ther- 
mal flux from the surface exerts a force on the current- 
carrying electrons in this region and a large electric 
field develops. A positron injected into this region 
quickly accelerates to energies where it can upscat- 
ter keV photons above the threshold for pair creation. 
The pairs emit synchrotron radiation consistent with 
the observed spectrum. 
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Figure 4: The left panel depicts the evolution of the effective temperature (left) and magnetic field (right) magnetic 
field decays through the solenoidal mode. The right panel depicts the results for decay through the irrotational mode. 
The bold solid, dotted and dashed lines give the results for Bi — 10^**, 10^^ and 10^® G respectively. The light lines give 
the results without any field decay, and the long-short-dashed line gives the cooling evolution for an unmagnetized 
neutron star. 



Heyl and Hernquist [70j present an alternative pic- 
ture that can account for both the optical and hard 
x-ray non-thermal emission. ii l3.1l outlines the dynam- 
ics of the electromagnetic field including radiative cor- 
rections. A key result is that the effective Lagrangian 
of QED makes electrodynamics non-linear. Specifi- 
cally simple electromagnetic waves travelling through 
a strong magnetic field develop regions when the elec- 
tric and magnetic field are discontinuous - shocks |^ . 
These shocks also form as fast waves travel through 
an ultramagnetized, ultrarelativistic plasma [t^ . 

As a fast wave propagates away from the surface of 
the neutron star and forms a shock, the energy release 
in the shock powers the formation of electron-positron 
pairs at rest in the frame of the shock. These pairs 
emit synchrotron photons that produce more pairs. 
This pair cascade stops because the total optical depth 
for pair creation strongly decreasing with the energy 
of the photons. The final generation of pairs has a 
typical Lorentz factor of 7 = 0.05^^^ with ^ ^ 1 
far from the surface of the star. These pairs emit 
synchrotron photons with an energy ranging from 

inbreak = J^nujB » 2.5 X IQ-^C^mc'. (33) 

down to Eq — ^mc^. At each radius the pairs cool 
in a fixed magnetic field, yielding a spectrum dE cx 

— 1/2 

E-y dEj between the values of Eq and £'brcak at the 
innermost edge of the pair production region. At en- 
ergies greater than inbreak or less than Eq the emission 
comes from pairs produced further from the star. Sim- 
ulations of the fast-mode breakdown in a dipole ge- 



ometry |70j indicate that the fast mode delivers equal 
amounts of energy in equal ranges of mangetic field, 
yielding the complete spectrum 

dE i ^min < E.y < Eq 

OC I Ey Eq < E^ < ii'bioak (34) 

[ -^7 ^-broak < Ej < E^^x 

The model has two free parameters: the total energy 
delivered by the fast mode i?totai and the strength 
of the magnetic field at the inner edge of the pair- 
production region that determines the location of both 
break points. 

Fig. [S] compares the spectrum of non-thermal 
radiation produce by fast-mode breakdown with 
the non-thermal radiation observed from the AXPs 
4U 0142+61 and IE 1841-045 and the SGR 1806-20. 
The model accounts for the observed spectral slope 
both in the optical and in the hard x-ray. If we as- 
sume that all of the objects have similar intrinsic op- 
tical emission to the nearby AXP 4U 0142-1-61, we 
predict that the non-thermal emission should peak 
around 100 MeV as traced by the curves denoted as 
the "Optical Model" for AXP 4U 0142-^61 and as the 
"Unified Model" for the other two objects. On the 
other hand if IE 1841-045 and the SGR 1806-20 do 
not exhibit a similar optical excess to 4U 0142-1-61, 
the non-thermal emission should extend to at least 
1 MeV as depicted by the "Minimal Model" curves. 

EGRET determined upper limits for the gamma- 
ray flux from the direction of these objects. Tab. |l] 
gives these upper limits and the predictions of the 
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Figure 5: The spectrum produced by fast-mode 
breakdown is superimposed over the observed thermal 
and non-thermal emission from several AXPs and SGRs 
for models that fit either the optical or INTEGRAL data 
solely and one that fits both sets of data. The 
unabsorbed optical data are from HuUeman et al. 
via Ozel 68] for AXP 4U 0142-1-61. The uppermost black 
symbols are the hard X-ray band are from Molkov et al. 
[S3 for SGR 1806-20. Mereghetti et al. ^ obtained 
similar results for the SGR. The middle sets of points in 
the hard X-ray data (blue is total flux and red is pulsed 
fiux) are from Kuiper et al.|2i| for AXP IE 1841-045. 
The green squares plot the INTEGRAL data reported by 
den Hartog et al. 0] for AXP 4U 0142~f61. We 
normalised the den Hartog et al.l6a| results using the 
observations of the Crab by Jung |75| . We scaled the 
emission from the three sources by assuming that they all 
lie at the distance of AXP 4U 0142-1-61. We used 3 kpc 
for AXP 4U 01424-61 7.5 kpc for AXP IE 1841-045 
Ili and 15 kpc for SGR 1806-20 113. 



fast-mode breakdown model. The optical model 
and the unified model predict approximately simi- 
lar EGRET fluxes for 4U 0142+61. The minimal 
model based solely on the INTEGRAL data exhibits a 
flux above 100 MeV about two hundred times smaller 
than the unified model. Because the optical model 
cannot explain the observed INTEGRAL data for 
IE 1841-045 and SGR 1806-20, it is omitted. Sim- 
ilarly, the minimal model cannot explain the opti- 
cal data for 4U 0142-f61. We see that for IE 1841- 
045 and SGR 1806-20 the predictions for the minimal 
model lie comfortably below the EGRET upper lim- 
its. In the context of the fast-mode breakdown model, 
this means that the optical emission for IE 1841- 
045 and SGR 1806-20 is inherently weaker than from 
4U 0142+61. 



On the other hand, 4U 0142+61 is difficult to ex- 
plain in the context of either model because of its 
large optical flux. Perhaps 4U 0142+61 was more ac- 
tive during the epoch of the optical observations than 
during the EGRET observations. Some AXPs exhibit 
variable X-ray emission such as AX J1845-0258 and 
IE 1048.1-5937 [11 [t^ so this conclusion might be 
natural. 



6. OUTLOOK 

Our understanding the magnetars has increased 
dramatically over the past few years but so have the 
unknowns. The recently discovered hard x-ray emis- 
sion ,6^ ,6^ |63, [3|P^^sent a theoretical challenge to 
understand but may also hold that key to understand- 
ing the evolution of the mangetic field on magnetars 
that drives not only the non-thermal and thermal qui- 
escent emission but the bursts as well. The recent 
massive burst from SGR 1806-20 [H IH HI shows 
that the size of soft-gamma repeater bursts varies over 
two additional orders of magnitude. The gamma-ray 
energy from the December 27 event (2 x lO''^^ erg) is 
a whopping two-percent of the total electromagnetic 
energy from a supernova, and the SGR is expected to 
burst like this many times over its lifetime. The ob- 
servation of SGR-like bursts from AXPs has further 
unified these two classes of objects and the association 
of glitches with the burstsf5i from AXP IE 2259+586 
gives a tantalizing hint at the underlying mechanism 
for these bursts. 

We still do not have a adequate theoretical descrip- 
tion of the thermal spectra from AXPs and SGRs. 
Where do the expected line features go? Alth ough 
features have been seen in IRXS J170849-400910 ifll, 
the other AXPs and isolated neutron stars such as 
RX J185635-3754 and RX J0720.4-3125 exhibit no fea- 
tures. Is this an indication that magnetars have con- 
densed surfaces 's^ or do they have helium, carbon 
or iron atmospheres that have not yet been studied in 
the magnetar regime? 

Understanding magnetars draws upon a wide range 
of physical processes from quantumelectrodynamics 
to plasma physics, from condensed matter physics 
to general relativity. The physics is messy, but that 
makes it fun! 
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Table I Predicted flux above 100 MeV and observed EGRET upper limi ts Fl in units of 10"* photons s"^cm"^ The 
GLAST upper limits are nominally (0.2 - 0.4) x 10"* photons s^^cm"^ l79l. 



Object 



EGRET EGRET Unified Minimal Optical 

Exposure [weeks] Upper Limit Model Model Model 



AXP 4U 0142+61 
AXP IE 1841-045 
SGR 1806-20 



8.8 
6.8 
4.9 



50 
70 
70 



1500 
70 

280 



800 



0.4 
0.6 
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